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Direct observation of the quantum critical point in heavy fermion CeRhSis 
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We report on muon spin rotation (/xSR) studies of the noncentrosymmetric heavy fermion an- 
tiferromagnet CeRhSis. A drastic and monotonia suppression of the internal fields, at the lowest 
measured temperature, was observed upon increase of external pressure. Our data suggest that 
the ordered moments are gradually quenched with increasing pressure, in a manner different from 
the pressure dependence of the Neel temperature. At 23.6 kbar the ordered magnetic moments are 
fully suppressed via a second order phase transition and Tn is zero. Thus, we directly observed the 
quantum critical point at 23.6 kbar hidden inside the superconducting phase of CeRhSis. 



PACS numbers: 76.75.+i, 75.50.Ee, 75.30.Mb, 75.40.Cx 



The role of magnetism in many superconducting ma- 
terials is a topic of intense research. Often, these mate- 
rials display a quantum critical point (QCP) that sepa- 
rates a nonmagnetic from a magnetic metallic phase fil- 
ial. This usually leads to the formation of new phases 
of matter, whose properties are dominated by quantum 
fluctuations. These fluctuations can be fine-tuned by an 
external parameter, such as chemical substitution, mag- 
netic field or external pressure. The location and nature 
of QCPs are central to the understanding of the inter- 
play between superconductivity and magnetism in these 
systems. In a number of heavy fermion (HF) compounds 
such as CePd2Si2, Celns, and CelrSia, the QCP is buried 
deep inside of a superconducting phase in the pressure- 
temperature phase diagram (1, 01 • In these cases it 
is challenging to establish the presence of a QCP and 
to study the microscopic magnetic properties coexisting 
with superconductivity. In this /iSR study on the HF 
CeRhSis, we uncovered the previously conjectured QCP 
hidden inside the superconducting phase. We directly 
observed the disappearance of the ordered magnetic mo- 
ments at the critical pressure p* = 23.6 kbar deep inside 
the superconducting phase. This suppression of the or- 
dered moments was found to be accompanied by a zeroing 
of the Neel temperature at p* suggesting that this QCP 
is of magnetic nature. 

CeRhSis is a noncentrosymmetric HF system jsl-lsj, 
suspected to display a QCP (sl, [ol, [lo[ at pressures of 
the order of 26 kbar. The compound features an intrigu- 
ing phase diagram. It displays an AFM phase at am- 
bient and applied pressure and a superconducting phase 
for pressures p > 12 kbar 0, Measurements of the 
de Haas- van Alphen effect at ambient and applied pres- 
sure show that the f -electrons are itinerant up to at 



least 30 kbar ^ . Neutron diffraction at ambient pres- 
sure showed that the magnetic order is incommensurate, 
with ordered moments of the order of 0.1 /ie/Ce (l2| . 
From previous results of thermal and transport proper- 
ties it was established that the antiferromagnetic order- 
ing temperature Tn ^ 1.6 K at ambient pressure, ini- 
tially increases with increasing pressure, and then de- 
creases towards the superconducting transition tempera- 
ture Tc ~ 1 K for p ^ 20 kbar. The nature of the mag- 
netism above this pressure across the proposed QCP in- 
side the superconducting phase is unclear. 

Single crystal CeRhSis platelets were grown via a stan- 
dard metal flux technique [l3| with initial molar ratios of 
1:1:3:30 (Ce:Rh:Si:Sn). The high quality of our single 
crystal samples was verified by the results of X-ray Laue 
diffraction and by measurements of the residual resis- 
tance ratio psoo k/po which was found to be of the order 
of 110. 

/iSR measurements at ambient pressure were per- 
formed using the Low Temperature Facility spectrometer 
(LTF) at the Paul Scherrer Institut (PSI). Measurements 
under pressure were conducted on the General Purpose 
Decay-Channel spectrometer (GPD). For zero-field (ZF) 
and longitudinal-field (LF) measurements we typically 
counted 6x 10^ positron events per spectrum and 0.5 x 10^ 
for transverse- field (TF) measurements, respectively. 

The single crystal sample consisted of platelets with 
an average thickness of 0.1 mm. For the measurement on 
LTF they were attached to a silver plate, with the crys- 
tal c-axis perpendicular to the plate and parallel to the 
muon beam covering an area of ~ 1 cm^. For pressure- 
measurements they were stacked in a cylindrical volume 
of 1 cm^ (mass ~ 500 mg). In addition, a 1 g polycrys- 
talline sample of CeRhSis was measured under pressure. 



2 




t{\LS) e (Gauss) 

FIG. 1. (Color online) ZF measurements at various ap- 
plied pressures on LTF (p = kbar) and GPD (p = 
6.8 and 16.9 kbar) from single crystal and polycrystalline 
CeRhSis samples, (a) The time spectra at different pressures 
at base temperature (0.02 and 0.27 K, respectively) and their 
corresponding fits (black lines) according to Eq. ([1)). (b) FFTs 
of the time spectra shown in (a). Both the spectra in (a) and 
the FFTs in (b) are vertically shifted for better illustration. 

In Fig. [1] ZF time spectra and the corresponding fast 
Fourier transforms FFTs for several pressures at a base 
temperature (of 0.02 K and 0.27 K for LTF and GPD, 
respectively), are shown for single crystal and polycrys- 
talline samples. The signal from the Ag plate and the 
pressure cell, respectively, has been subtracted, (see Sup- 
plemental Material for details). The signals shown in 
Fig. [U reveal a fast decay and a characteristic time struc- 
ture, expected in the case of a signal containing sev- 
eral oscillatory components. The fast Fourier transform 
(FFT) of the time spectrum at p = kbar (see Fig.[T](b)), 
displays three well-defined distributions of internal mag- 
netic fields centered at ~ 0, 45 and 120 G. Therefore, 
the ZF sample asymmetry function is described here by 
three terms: 
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gi{t) = a • cos (j^B.t) • e"^^^ ' ^ + (1 - a) • e'^^ '\ 

where 7/i/27r = 135.5 MHz/T with 7^ the muon gyro- 
magnetic ratio. Bi denote the average internal magnetic 
fields at the muon stopping sites in the sample, a was 
found to be | , which is not apriori evident for single crys- 
tals, but is always the case for polycrystalline samples. 
Ai are the asymmetries at time of the individual com- 
ponents, which are proportional to the probability that 
a muon stops at a site with an internal field Bi. We find: 
Ai = A2 and As = 1.5 • Ai. Axi and Al are the transverse 
and longitudinal muon relaxation rates, respectively. Bi, 
Axi and Al are temperature dependent. The relations 
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FIG. 2. (Color online) ZF data from the polycrystalline sam- 
ple at (a) 20.4 kbar and (b) 23.6 kbar at two different tem- 
peratures. The time spectra were fitted using Eq. ([Tj). The 
signal of the pressure cell was subtracted. 



B1/B2 ^ 2.6 and Axi/Ax2 ^ 2.6 were found experimen- 
tally; Bs was found to be roughly for all temperatures 
and pressures. 

The weights of the spectra in Fig. [TJb) shift to lower 
fields with increasing pressure, that is, the quasi-static 
internal fields are increasingly suppressed with pressure. 

Since the AFM structure is incommensurate, one 
might expect that the quality of the fits would improve 
by replacing the cosine terms in Eq. ffl^ by Bessel func- 
tions, as in the case of UNi2Al3 [ij, E|, CeAls ^ or 
CeCusAu This replacement, however, leads to fits 
of lower quality, and will not be considered further here. 

High-pressure data at low and high temperatures 
are shown in Fig. [21 In the time spectrum at 20.4 kbar 
and T = 0.27K <C Tn we observe a fast relaxation and 
slow oscillations. Both of these features are absent at 
T = 1.76 K > Tn showing that the low-temperature data 
reflects the presence of internal magnetic fields arising 
from magnetic order. For the data at 23.6 kbar, the fast 
relaxation and the oscillations are absent and there is no 
significant change from the high temperature to the low 
temperature signal at T = 0.27 K. This demonstrates 
the absence of magnetic order down to at least 0.27 K, 
the lowest temperature of these experiments. Since at 
this pressure the internal magnetic fields are zero, Tn is 
most likely zero. 

In the results of resistivity and specific heat measure- 
ments CeRhSis was found to be superconducting for 
p > 12 kbar and low T, our results therefore demonstrate 
that over the pressure region 12 kbar < p < 23.6 kbar su- 
perconductivity and AFM order coexist, as indicated in 
earlier measurements jij. 

The ZF data in the AFM state shows, in addition to 
the tiny ordered moments, a fast initial relaxation. As 
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FIG. 3. (Color online) /iSR time spectra of the single crys- 
tal sample in zero and longitudinal field at T — 0.27 K and 
p = 16.3 kbar. The pressure cell signal (background) was 
subtracted. 

we shall see, the fast relaxation is not caused by fluctu- 
ating (large) moments, as it would be in the case where 
Kondo screening were destroyed and the local Ce mo- 
ments were restored. The results of /iSR measurements 
at longitudinal fields at p = 16.3 kbar (see Fig. [3]) showed 
that at modest external fields of 200 gauss (of the order 
of the internal quasi-static fields) the muon relaxation 
is substantially suppressed; that is, magnetism is fully 
decoupled, which supports the quasi static origin of the 
relaxation at the pressure 16.3 kbar. We find no evi- 
dence for dynamical relaxation at any pressure in the 
low-temperature region of the AFM phase. Above Tn at 
p = kbar an LF measurement was also conducted and 
complete decoupling was observed at an applied field of 
~ 200 G. We therefore infer absence of appreciable dy- 
namical relaxation above Tn- Thus, our data supports a 
scenario where the Ce moments are quenched at all pres- 
sures below p*, and where the Kondo screening remains 
largely intact. 

The normalized static internal fields Bi/ Bi{0){p) 
(equivalent to B2/ B2{0){p)) extracted from the fits using 
Eq. ([1]) are depicted in Fig. [H as a function of pressure 
(circles). Bi{0) denotes the internal field at the lowest 
measured pressure p = kbar. Bi decreases monoton- 
ically with pressure and vanishes near p* ~ 23.6 kbar, 
showing that magnetic order is suppressed at . For the 
two highest measured pressures the field distribution is 
very broad. Still, over the entire antiferromagnetically 
ordered region one observes a pressure independent ratio 
B1/B2. This implies that the magnetic structure does 
not change appreciably with pressure. In addition, our 
data show that near p* the fast relaxation characterizing 
the ZF data in the AFM state vanishes. 

The Neel temperature as a function of pressure was 



p (kbar) 

FIG. 4. (Color online) Comparison of the pressure depen- 
dence of the Neel temperature (diamonds) and the normal- 
ized internal field Bi / Bi (0) (circles) . The lines are guides to 
the eye. The inset shows the asymmetry Ai^s versus temper- 
ature exemplarily for the pressure p = 20 A kbar. The Neel 
temperature is indicated by an arrow. The dashed line is a 
guide to the eye. 

obtained with /iSR measurements in transverse fields of 
50 and 300 gauss, respectively. Below the Neel tempera- 
ture the magnetic order leads to an extra fast relaxation 
so that the total asymmetry after some time, t > Ifis 
in our case, is smaller than the corresponding amplitude 
above Tn- The time spectra for t > 1 /is were fitted using 

A{t) = Ai^s • exp ( - ^icTtf) • cos {-/^Bt), (2) 

where Ai^g is the asymmetry at 1 /is after the muon 
is implanted. B represents the local field and is a fit 
parameter. Upon cooling, one observes a rapid decrease 
of as the system enters the AFM state. The 

Neel temperature is associated with the onset of the 
drastic steplike reduction of Ai^s{T) with decreasing 
temperature. One example of Ai^s{T) for p = 20.4 kbar 
is shown in the inset of Fig. IH 

The results for Tn at different pressures are shown in 
Fig. [Has diamonds. Tn(p) shows a maximum at around 
7 kbar, consistent with the results of resistivity (sl, [llj 
and specific heat measurements. By contrast, the 
internal fields and, therefore, the ordered magnetic 
moments decrease monotonically. We find that the 
internal fields and Tn vanish at the critical pressure of 
p* = 23.6 kbar and magnetic order disappears. Our 
data show that the AFM order is completely suppressed 
for p > p*, and that there is a continuous tuning of 
the AFM phase with pressure. All this indicates the 
presence of a magnetic QCP at p*. 

The different pressure dependencies of the Neel 
temperature and that of the internal fields may be 
understood using the simplified, but widely used. 
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one-dimensional Kondo lattice model described by S. 
Doniach fll|. The model considers the competing effects 
of the Kondo screening and the exchange interaction 
between Ce moments which tend to stabilize the mag- 
netic moments leading to magnetic order. The results 
are summarized in the so-called Doniach phase diagram, 
displaying the Neel temperature Tn as a function of 
the product |JA^f| of the exchange coupling J and the 
density of states of the conduction electrons at the 
Fermi level N^. For low values of |JA^f| the system is 
magnetically ordered and Tn increases with increasing 
|JA^f|- Further increase of |JA^f| yields a maximum of 
Tn and a subsequent decrease (observed, for instance 
in CeAg Tn eventually reaches zero, and further 

increase in |JA^f| yields a nonmagnetic (Fermi liquid) 
system. 

In our case, the product \JNy\ is tuned with external 
pressure. The observation of the maximum of T^{p) 
and the subsequent suppression of Tn to zero, suggests 
that the application of rather modest pressures leads 
to substantial changes of \JNy\. At the same time, the 
Doniach model predicts a monotonic decrease of the 
ordered moments with increasing J according to [3l 
{S^) = y/l — (J/W)2, where W is the conduction 
electron bandwidth and is the x-component of a 
spin S. These predictions match qualitatively our 
experimental observations. 

We emphasize that the order parameter, i.e., the sub- 
lattice magnetization (or more precisely the amplitude of 
the spin density wave), represented here by the internal 
fields, displays a continuous phase transition with 
pressure and the Neel temperature goes to zero at p*. 
This clearly reveals the presence of a QCP. In our case, 
it lies deep in the superconducting phase, and its bulk 
properties are masked by superconductivity. Neverthe- 
less, it has been reported that between 20 — 26 kbar and 
T > 1 K, no n- Fermi-liquid behavior [6] is observed in 
the temperature dependence of the electrical resistivity, 
adding support to our claim. Further indications for a 
QCP may be extracted from the extremely high upper 
critical field Hc2 (Hc2\\c) and its sensitivity to pressure 
for p > 26 kbar [20| . Both the behavior of the resistivity 
and that of Hc2 may be explained theoretically in 
terms of antiferromagnetic spin fluctuations in three 
dimensions close to the QCP [lo|- In the results of de 
Haas-van Alphen measurements no change in the size 
of the Fermi surface across p* was observed [9]. That 
is, the Kondo screening of the Ce moments is active at 
pressures on both sides of p*. All this hints at a SDW 
type QCP [H, consistent with our findings. This 
type of QCP is driven by quantum fluctuations of the 
magnetic order parameter, as opposed to the so called 
local QCP where the Kondo screening is destroyed and 
the local moments are restored in the magnetically 
ordered phase, i.e. p < p^. 

In summary, we have measured the pressure evolution 



of the internal fields of the AFM phase of CeRhSis 
with the /iSR technique. As a function of pressure the 
internal fields diminish continuously with increasing 
pressure, in clear contrast to the complex behavior of the 
Neel temperature which shows a maximum at around 
7 kbar. Our observations may be understood in the 
framework of the Doniach phase diagram. This analysis 
reveals a strong sensitivity of \JNy\ to the externally 
applied pressure for CeRhSis. At p* ^ 23.6 kbar the 
system becomes non-magnetic and signals the presence 
of a magnetic QCP, masked by the superconducting 
phase. We find no evidence for a breakdown of the 
Kondo screening across the QCP. 
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Supplemental material: Background signal 

In our experiments the raw /iSR signal contains two com- 
ponents: the signal from the sample, CeRhSis, and the 
signal from the background (bg), i.e. 

A(t) = AceRhSi3(^)+^bg(t). (3) 

In LTF the background signal originates from muons 
stopped in the silver plate to which the crystals are 
attached, whereas in GPD the background stems from 
muons stopped in the pressure cell. Both cases have 
been extensively characterized in previous experiments 
and their relaxation is due to the nuclear moments of the 
materials. In the analysis of our data, the relaxation rates 
of the respective backgrounds were fixed in the fits and 
the only free parameter is the initial asymmetry A\)g{0). 
With this restriction we perform a full fit to the data us- 
ing Eq. m from where we extract AceRhSis • 

The silver plate used in the LTF experiment con- 
tributes a signal that is roughly constant in the measured 
time window. Approximately 50% of muons stopped in- 
side of the sample in the LTF measurements. 

For the GPD measurements, we used piston-cylinder 
pressure cells made of CuBe {p < 21 kbar) and the Co- 
Ni alloy MP35 {p > 21 kbar), respectively (2l|. The hy- 
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drostatic pressure was measured with an 0.1 kbar accu- 
racy in the fohowing way: with ac-susceptibihty measure- 
ments the superconducting transition temperature of 
an indium piece, placed inside the pressure ceh, was de- 
termined. The pressure dependence of the superconduct- 
ing transition of In is 



Tc = 3.402 K- 3.64.10"^ K/kbar-p, 



(4) 



where p is the pressure in kbars [22| . 

The pressure ceh contributes a rather high background 
compared to the signal from the sample. For the sin- 
gle crystal sample only 20 - 25% of the incoming muons 
stopped in the sample. For the poly crystalline sample 
the stopping fraction increased to 40%. 

In zero field (ZF) the signal due to the CuBe-pressure 
cell can be described by a Kubo-Toyabe type polarization 
function [Is^: 



^ktW = - + -(1 



(at)2)exp[--(at)2], (5) 



where is the width of the local field distribution 

with <j = 0.33 /is~^ (temperature independent) for CuBe. 

For the MP35-pressure cell the signal in ZF may be 
approximated by the product of a Kubo-Toyabe (Eq. (j5j)) 
and an exponential function, i.e. (2I, 24| 



^MP35(ZF)W = ^ktW • exp(-At). 



(6) 



The relaxation rate A is found to be temperature depen- 
dent and was fixed at each temperature to values from a 
calibration curve. The calibration curve was obtained at 
the beginning of the experiment at the measured pres- 
sure. For this purpose the energy of the incoming muons 
is adapted so that the muons stop only inside of the pres- 
sure cell (and not in the sample). In a weak transverse 
field (TF) the pressure cell is fitted by 

^MP35(TF) = ^0 exp[-i(a t)^] exp(-At) cos {-f^Bt), 

(7) 

with Aq the asymmetry at time zero, a' ^ 0.33 /is ^ 
(temperature independent) and B the local field parame- 
ter. A(T) is marked by a strong increase below T = 1.4 K. 
The subtraction of the MP35-pressure cell signal in the 
low temperature region therefore requires special care, 
due to its strong temperature dependence. 



The signal of the (CuBe) pressure cell in a longitudinal 
field (LF) is described by an extended Kubo-Toyabe type 
function [Eq. 21 of Ref. @ and Eq. 6.130 of Ref. ^]: 



^LF(t) = 1 - ^^[cOSC^t - jo(c^t)] 



■L 



t nil 



Q\s)fo{ujs) 



ds, 



(8) 



with jo{x) = sin{x)/x and Q{t) = exp(— ^a^t^). 
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